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\ We study conservation laws and potential symmetries of (systems of) differential equations 

applying equivalence relations generated by point transformations between the equations. 
A Fokker-Planck equation and the Burgers equation are considered as examples. Using 
reducibility of them to the one-dimensional linear heat equation, we construct complete 
hierarchies of local and potential conservation laws for them and describe, in some sense, all 
their potential symmetries. Known results on the subject are interpreted in the proposed 
framework. This paper is an extended comment on the paper of J.-q. Mei and H.-q. Zhang 
CO ■ [Internat. J. Theoret. Phys., 2006, in press]. 



^ ! 1 Introduction 

m . 

A number of authors investigate symmetries (in particular, potential ones) and conservation laws 
of different classes of equivalent equations. In such way they perform a huge number of unnec- 
essary cumbersome calculations. Since equivalence transformations may be quite complicated, 
sometimes it seems to be impossible to obtain directly complete and correct results for some 
cases which are equivalent to simple ones, although these results can be easily reconstructed 
Qh' with application of equivalence transformations. 

At the same time, it is a general mathematical rule that equivalent in some sense objects pos- 
sess equivalent properties. In particular, if two systems of equations are equivalent with respect 
to point transformations then there exists a one-to-one correspondence between their maximal 
Lie invariance algebras, spaces of conservation laws, potential symmetries, exact solutions etc. 
All the above mentioned features for more complicated models can be constructed from ones 
of the simpler equivalent models by means of application of known equivalence transformation 
■ easier than by direct calculations. 

In the presented paper we apply this observation to describe exhaustively local and potential 
conservation laws and potential symmetries of a Fokker-Planck equation and the famous Burgers 
equation due to reducibility of them to the linear heat equation. 

This paper was started as a comment on the recent work [14], where the authors tried to 
obtain some potential conservation laws of the Burgers and Fokker-Planck equations, referring 
to results of [19,20]. At first our aim was to correct mathematically and historically inaccurate 
statements of [14]. Thus, e.g., system (13) v x = u, vt = u x + xu is not a potential system for 
equation (12) u t = (xu x ) x + u xx of [14]. Comparing (12) and (13) with [19,20] we deduce that 
the authors consider conservation laws of the Fokker-Planck equation of form ut = u xx + (xu) x . 
Let us mention also that the potential operator (14) given in [14] with a reference to [19, 20] is 
written with missprints (up to the multiplier e 2< ). Therefore, the conservation law obtained with 
this operator may be incorrect. Potential symmetry operators for the Burgers equation obtained 
in [14] do not form linear space and, therefore, do not form a Lie algebra. Direct calculations 
show also that (T 1 ,!" 2 ) adduced in the last formula of the paper [14] is not a conserved vector 
for the Burgers equation. While working on the comment we decided to generalize essentially 
results of [14] adducing also the complete sets of the potential conservation laws and giving a 
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complete, in some sense, description of all possible potential symmetries of the Burgers and 
Fokker-Planck equations. 

To construct non-local conservation laws of the Fokker-Planck and Burgers equations the 
following statement [9,12] was used in [14] (see also [8,13,17]). 

Proposition 1. If the system C admits a one-parameter group of symmetry transformations with 
the infinitesimal generator Q = Tdt + £d x + r]d u and a conservation law of form D(T + D X X = 0, 
then the vector-function (T, X) with 

f = Q{T) + TD x i - XD x t, X = Q(X) + XD t r - TD£ 

is a conserved vector. 

If the conservation law is invariant with respect to the action of operator Q then (T, X) van- 
ishes. Thus, one obtains a rule for finding conservation laws invariant with respect to symmetry 
operator: 

Q(T) + TD x ti - XD x t = 0, Q(X) + XD t r - TD£ = 0. 

Although the latter formula (Theorem 3.1 from [14]) is completely true and can be useful for 
obtaining particular classes of conservation laws, it is not guarantee that the resulting conserva- 
tion law will be non-trivial or pure nonlocal. Thus, e.g., if one redo the calculations leading to 
formula (26) of [14] without computational inaccuracies, one obtains exactly conservation law 
equivalent to local one. In general, invariant conservation laws do not generate the complete 
space of conservation laws. Therefore, this approach is suitable only in case the complete set 
of conservation laws is difficult or even impossible to be constructed. As shown in [17], com- 
plete hierarchies of potential conservation laws for both equations under consideration can be 
constructed. 

Below we apply the equivalence framework to examine potential symmetries and conservation 
laws of the Fokker-Planck and Burgers equations and to give a detailed interpretation of inter- 
esting results of Pucci and Saccomandi [19,20]. Since the both equations under consideration 
are reduced to the linear heat equation, we also describe a way to find potential symmetries of 
the linear heat equation using potential systems associated with nonconstant characteristics. 

More precisely, our paper is organized as follows. At first, in Section [2] we adduce a necessary 
theoretical background on conservation laws and potential systems, including a short discussion 
of empiric and rigorous definitions of conservation laws. Equivalences of conservation laws and 
potential systems with respect to point transformations are also discussed in detail. Then, using 
reducibility of the Fokker-Planck and Burgers equations to the one-dimensional linear heat 
equation, we construct complete hierarchies of local and potential conservation laws for them 
and describe their potential symmetries (Sections [3] and [5]). Potential symmetries of the linear 
heat equation are considered in Section HI 

2 Theoretical background 

To begin with, we adduce a necessary theoretical background on conservation laws and potential 
symmetries, following, e.g., [16,17,24] and considering for simplicity the case of two independent 
variables t (the time variable) and x (the space variable). See the above references for the general 
case. 

Let C be a system L(t, x, u/ p )) = of / PDEs L l = 0, . . . , L l = for the unknown functions 
u = (u\ . . . , u m ) of the independent variables t and x. Here ur p \ denotes the set of all the partial 
derivatives of the functions u of order not greater than p, including u as the derivatives of the 
zero order. 

First we give an empiric definition of conservation laws. 
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Definition 1. A conservation law of the system £ is a divergence expression 



D t T(t,x,U( r )) + D x X(t,x,u^) = 



(1) 



which vanishes for all solutions of £. Here Dt and D x are the operators of total differentiation 
with respect to t and x, respectively. The differential functions T and X are correspondingly 
called a density and a flux of the conservation law and the tuple (T, X) is a conserved vector of 
the conservation law. 

The crucial notion of the theory of conservation laws is one of equivalence and triviality of 
conservation laws. 

Definition 2. Two conserved vectors (T, X) and (T", X') are equivalent if there exist func- 
tions T, X and H of t, x and derivatives of ?/ such that T and X vanish for all solutions of £ and 
T' = T + T + D X H, X' = X + X — D t H. A conserved vector is called trivial if it is equivalent 
to the zeroth vector. 

The notion of linear dependence of conserved vectors is introduced in a similar way. Namely, a 
set of conserved vectors is linearly dependent if a linear combination of them is a trivial conserved 
vector. 

It is obvious that under the problem of finding conservation laws for some system one should 
understand the problem of finding inequivalent linearly independent conservation laws, i.e., con- 
servation laws having linearly independent conserved vectors. 

Conservation laws can be investigated in the above empiric framework. However, for deeper 
understanding of the problem and absolutely correct calculations a more rigorous definition of 
conservation laws should be used. 

For any system £ of differential equations the set CV(£) of conserved vectors of its con- 
servation laws is a linear space, and the subset CVo(£) of trivial conserved vectors is a linear 
subspace in CV(£). The factor space CL(£) = CV(£)/ CVo(£) coincides with the set of equiv- 
alence classes of CV(£) with respect to the equivalence relation adduced in Definition [2j 

Definition 3. The elements of CL(£) are called conservation laws of the system £, and the 
whole factor space CL(£) is called the space of conservation laws of £. 

That is why description of the set of conservation laws can be assumed as finding CL(£) 
that is equivalent to construction of either a basis if dimCL(£) < oo or a system of generatrices 
in the infinite dimensional case. The elements of CV(£) which belong to the same equivalence 
class giving a conservation law T are considered all as conserved vectors of this conservation law, 
and we will additionally identify elements from CL(£) with their representatives in CV(£). For 
(T, X) G CV(£) and T € CL(£) the notation (T,X) € T will denote that (T,X) is a conserved 
vector corresponding to the conservation law T . In contrast to the order rrp x) °f a conserved 
vector (T, X) as the maximal order of derivatives explicitly appearing in the differential functions 
T and X, the order of the conservation law T is called min{rrr,x) I (T,X) G J 7 }. Under linear 
dependence of conservation laws we understand linear dependence of them as elements of CL(£). 
Therefore, in the framework of "representative" approach conservation laws of a system C are 
considered as linearly dependent if there exists linear combination of their representatives, which 
is a trivial conserved vector. 

Let the system £ be totally nondegenerate [16]. Then application of the Hadamard lemma to 
the definition of conservation law and integrating by parts imply that the left hand side of any 
conservation law of £ can be always presented up to the equivalence relation as a linear combi- 
nation of left hand sides of independent equations from £ with coefficients being functions 
of t, x and derivatives of u: 



D t T + D X X = ^L 1 + ■■■ + \ l L l . 



(2) 
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Definition 4. Formula ([2]) and the £-tuple A = (A 1 , . . . , A') are called the characteristic form 
and the characteristic of the conservation law D{T + D X X = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of C. Since L is nondegenerate, 
the characteristics A and A satisfy © for the same conserved vector (T, X) and, therefore, 
are called equivalent iff A — A is a trivial characteristic. Similarly to conserved vectors, the 
set Ch(£) of characteristics corresponding to conservation laws of the system C is a linear 
space, and the subset Cho(£) of trivial characteristics is a linear subspace in Ch(£). The factor 
space Chf(£) = Ch(£)/ Cho(£) coincides with the set of equivalence classes of Ch(£) with 
respect to the above characteristic equivalence relation. 

Any conservation law ([T]) of jC allows us to deduce the new dependent (potential) variable v 
by means of the equations 

v x = T, v t = -X. (3) 

In the case of single equation C, equations of form ([3|) combine into the complete potential 
system since C is a differential consequence of ([3]). As a rule, systems of such kind admit a 
number of nontrivial symmetries and so they are of a great interest. If the transformation 
of some of nonlocal variables t, x or u depends explicitly on variable v, such symmetry is a 
nonlocal for the initial equation (system) and is called potential symmetry. Let us mention that 
the concept of potential symmetry was introduced by Bluman at al [5, 6] in the late 80-es. See 
also the related notion of quasilocal symmetry [1,2]. 

An important property of the class of equations in the conserved form is that it is preserved 
under any point transformation (see, e.g., [17]). 

Proposition 2. A point transformation g: t = t 9 (t,x,u), x = x 9 (t,x,u), u = u 9 (t,x,u) 
prolonged to derivatives of u transforms the equation D{T + D X X = to the equation 
DfT 9 + D X X 9 = 0. The transformed conserved vector (T 9 , X 9 ) is determined by the formula 

. . T(x,ur r ))D t i+X(x,v,M)D x i 

T 9 (x,U( r \) = L 4 = — — , 

v ( >' D t t D x x - D x t D t x 

_ T(x,ur r ))D t x + X(x,ur r ))D x x 

X 9 (x,U( r )) = = — — . 

v (>J D t t D x x - D x t D t x 

Note 1. In the case of one dependent variable (m = 1) g can be a contact transformation: 
t = t 9 (t, x, U(i)), x = x 9 (t, x, M(i)), = u 9 ^(t, x,U(i^). Similar notes are also true for the 
below statements. 

Proposition 3. Any point transformation g between systems C and C induces a linear one- 
to-one mapping g* from CV(£) into CV(C), which maps CYq(C) into CVo(£) and generates a 
linear one-to-one mapping g^ from CL(£) into CL(£). 

Corollary 1. Any point transformation g between systems C and C induces a linear one-to-one 
mapping g^ from Chf(£) into Chf(£). 

It is possible to obtain an explicit formula for correspondence between characteristics of C 
and C. Let & = A^ U C U , where = A» ua D a , A^ ua are differential functions, a = (a t ,a x ) 
runs the multi-indices set (a^ GN U {0}), [A, v = 1, 1. Then 

A M = A^*((D t tD x x - D x tD t x)X u ). 

Here A^* = (— D) a ■ A tiVa is the adjoint to the operator A Vfi . For a number of cases, e.g. if C 
and C are single partial differential equations (1 = 1), the operators A^ v are simply differential 
functions (i.e., A^ ua = for \a\ > 0) and, therefore, A^* = A^ u . 

Equivalent conservation laws give rise to equivalent potential systems. More precisely, Propo- 
sition [2] and equations © imply the following statement. 
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Proposition 4. Any point transformation connecting two systems C and C of PDEs with two 
independent variables generates a one-to-one mapping between the sets of potential systems, 
which correspond to C and C Generation is made via trivial prolongation on the space of 
introduced potential variables, i.e., we can assume that the potentials are not transformed. 

In such way, if a transformation connects two systems of differential equations, then the same 
transformation maps the set of conservation laws of the first system to the set of conservation 
laws of the second system and the space of characteristics of conservation laws to the space 
of characteristics. This transformation is trivially prolonged to the possible potential variables 
and then makes a mapping between the corresponding sets of potential systems. Therefore, one 
can easily derive characteristics of conservation laws, conservation laws, potential systems and 
potential symmetries of the more complicated system from the ones of the simpler system. 



3 Fokker— Planck equation 

In [19,20] Pucci and Saccomandi investigated potential symmetries of the Fokker-Planck equa- 
tion 

u t = u xx + (xu) x . (4) 

Using the simplest conservation law DfU — D x (u x + xu) = with the characteristic 1, they 
construct the corresponding potential system 

v x = u, v t = u x + xu (5) 

and then found its maximal Lie invariance algebra 

01 = (d t , e~ l d x , e' 2t d t - e~ 2t xd x + e~ 2t ud u , e l d x - e\xu + v)d u - e*xvd v , 

e 2t d t + e 2t xd x - e 2t (x 2 u + 2xv + 2u)d u - e 2t {x 2 + l)vd v , ud u + vd v , f x d u + fd v ), 

where the function / = f(t, x) runs the solution set of the equation ft = f xx + xf x . Any operator 
from gi is a potential symmetry operator of equation It is a nontrivial potential symmetry 
operator iff the coefficient of d u depends on v. Let us note that the algebra gi of potential 
symmetry operators differs from the maximal Lie invariance algebra 

g = (d t , e~ l d Xl e~ 2t d t - e~ 2t xd x + e~ 2t ud u , e t d x - e l xud u , 
e 2t d t + e 2t xd x - e 2t x 2 ud u , ud u , f x d u + fd v ) 

of equation (Jl|) and is not projectible to go- At the same time, the algebras go and gi are 
isomorphic. 

We interpret the results by Pucci and Saccomandi on the characteristic 1 in the above frame- 
work of the equivalence relation between conservation laws of different equations, which is ex- 
tended to equivalence of potential systems. Then we give a generalization for the case of arbitrary 
characteristic. 

It is well-known that the Fokker-Planck equation ((4|) is reduced to the linear heat equation 
u>i = u xx by the transformation 

T: t = -e 2t , x = e l x, u = e~ t u. 

The same transformation maps the conservation law of the Fokker-Planck equation with char- 
acteristic 1 to the one of the linear heat equation with the same characteristic. In view of 
Proposition U potential system ([5]) reduces to the potential system 

v x = u, vi = Ux (6) 
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of the linear heat equation by the transformation T trivially prolonged to the potential variable 
v = v. Therefore, the prolonged transformation T pT establishes an isomorphism between the 
algebra gi and maximal Lie invariance algebra 

01 = (dt, d x , 2td t + xd x - ud u , 2td x - (xu + v)d u - xvd v , 

it 2 dt + 4txd x — ((x 2 + 6t)u + 2xv)d u — (x 2 + 2t)vd v , ud u + vd v , h x d u + hd v ), 

of system ([6]) in a similar way as the simple transformation T do this between the algebra go 
and maximal Lie invariance algebra 

g = (d t , d x , 2td t + xd x , 2td x - xud u , At 2 d t + 4txd x - (x 2 + 2t)ud u , ud u , hd u ) 

of the linear heat equation. Here the function h = h(t, x) runs the solution set of this equation. 
Tildes over variables in the formulas for go and gi are omitted. 

As shown in [18], the algebra gi also is isomorphic to the maximal Lie invariance algebra of 
the potential equation Vf = v xx (the form of which coincides with the initial equation on u). 
Moreover, the Lie symmetry operators of the potential system © are the first prolongations of 
the Lie symmetry operators of the potential equation. It means that a similar statement is true 
for all equivalent equations. Namely, the Lie symmetry operators of the potential system ((^J) are 
the first prolongations of the Lie symmetry operators of the corresponding potential equation 
vt = Vxx + xv x that can be reduced to the "potential" linear heat equation = v xx by means 
of the truncated transformation t = \e 2t , x = e t x, v = v. 

The set of all possible linearly independent local conservation laws of the linear heat equation 
is well-known [7, 10, 11, 16, 17, 22] and consists of ones having the form 

D^au) + D x (a x u - au x ) = 0, 

where a = ot(t, x) is an arbitrary solution of the backward linear heat equation 6^ + a xx = 0. 
Using Proposition [2] or different versions of the direct method for finding conservation laws [3, 
4,17,23], we obtain the set of linearly independent local conservation laws of the Fokker-Planck 
equation (jlj) 

D t (au) + D x ((a x - xa)u - au x ) = 0, (7) 

where a = a(t, x) is an arbitrary solution of the linear equation at + a xx — xa x = 0, which is 
adjoint to the Fokker-Planck equation. 

In [17] the theorem was proved that any potential conservation law of the linear heat equation 
is equivalent to the local one. Therefore, in view of Proposition [2] the same is true for the Fokker- 
Planck equation (j3|), and formula ([7]) gives the complete description of the local and potential 
conservation laws of the Fokker-Planck equation (|4|) . 

Let us emphasize that the above statement on the potential conservation laws is true not only 
for the laws obtained from the potential system corresponding to the characteristic a = 1, but 
also for the system which is a union of any finite number of the potential systems corresponding 
to linearly independent solutions of the backward heat equation. 

4 On potential symmetries of linear heat equation 

In contrast to the conservation laws, potential symmetries of the linear heat equation are inves- 
tigated only in case of the single characteristics 1. The problem of construction of all possible 
potential symmetries of the linear heat equation remains open (in particular, it includes inves- 
tigation of systems with arbitrary families of linearly independent characteristics) . 
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Here we study the families of potential systems of the linear heat equation which are con- 
structed with single local conservation laws and make a preliminary investigation. As mentioned 
in Section [3l any local conservation law of the linear heat equation ut = u xx has the form 

D t (au) + D x (a x u - au x ) = 0, 

where a = a(t, x) is an arbitrary non-zero solution of the backward heat equation at + a xx = 0. 
The characteristic of the conservation law coincides with a. The associated potential system is 

v x = au, v t = au x — a x u. (8) 

The initial equation on u is a differential consequence of system ([8]). Another differential conse- 
quence is the equation 

v t + 2— v x -v xx = (9) 
a 

on the potential dependent variable v which is called the potential equation associated with the 
equation ut = u xx and the characteristic a. 

Consider a Lie symmetry operator Q = rd t + £d x + r\d u + 9d v of system ([8]). The coefficients 
of Q are functions of t, x, u and v. The infinitesimal invariance criterion [15, 16] implies for 
system (JSJ) that, in particular, 

T u = tu = u = 0, (10) 
r x = r v = = 9 VV = 0, (11) 

r?=k-e*--r-^eW^. (12) 

V a a / a 

The subsystem (jlOp of determining equations means that any Lie symmetry transformation 
of ([8]) with respect to t, x and v does not depend on u. Equation ([9]) is a differential conse- 
quence of system ([8]), and there exist one-to-one correspondence between the sets of solutions 
of equation ([9]) and system ([8]). Therefore, the truncated operator Q = rdt + £,d x + 6d v is a Lie 
symmetry operator of equation ([9]). 

And vice versa, consider a Lie symmetry operator Q = rdt + £,d x + 9d v of equation ([9]). The 
coefficients of Q are functions of t, x and v. Then the prolonged to u operator Q = Q + r]d u , 
where r/ is defined by formula (|12p . is a Lie symmetry operator of system (jHJ). 

In view of the subsystem (jlip of determining equations and the formula for rj, we have 
Vv = 6 xv Hence the conservation law with the characteristic a results to pure potential symme- 
tries of the linear heat equation iff the potential equation possesses a Lie symmetry operator 
with the coefficient 9 of d v , which satisfy the condition 9 XV ^ 0. 

Let us study the case of the simplest nonconstant characteristic a = x. The corresponding 
potential equation has the form 

2 . . 

v t + -v x - v xx = 0. (13) 

x 

The maximal Lie invariance algebra of (|13|) is 

t = (d t , 2td t + xd x , At 2 d t + Atxd x - (x 2 - 2t)vd v , vd v , fd v ), 

where the function / = f(t,x) runs the solution set of (|13p . It is easy to see that for the third 
basis operator the condition 9 XV ^ is satisfied. Therefore, consideration with the characteristic 
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a = x leads to pure potential symmetries of the linear heat equation. Namely, the system (|8|) 
with a = x possesses the maximal Lie invariance algebra 

ti = (d t , 2td t + xd x - 2ud u , 4t 2 d t + 4txd x - {{x 2 + 6t)u + 2v)d u - {x 2 - 2t)vd v , 
ud u + vd v , x~ l f x d u + fd v ), 

where the function / = f(t,x) again runs the solution set of (|13p . Any linear combination 
of operators from ti which contains the third basis operator is an pure potential symmetry 
operator of the linear heat equation. Note that the pure potential symmetry operators from 
the algebra ti differ from ones from the algebra Qi both the explicit form and the nature of the 
potential variable which is defined by system (|8|) with a = x instead of system ([6]). 

5 Burgers equation 

The famous Burgers equation 

ut = u xx + 2uu x (14) 

is singular among nonlinear evolution equations due to its symmetry properties and admits the 
five-dimensional maximal Lie invariance algebra [16] 

b = (d t , d x , td x - d u , 2td t + xd x - ud u , t 2 d t + txd x - (tu + x)d u ). 

It has one linearly independent local conservation law DfU — D x (u x + u 2 ) = 0, the characteristic 
of which equals to 1. The Lie symmetries of the potential system 

v x = u, v t = u x + u 2 

associated with this conservation law are also well-known and studied by many authors. See, 
e.g., [18,21]. Its maximal Lie invariance algebra 

t>i = (dt, d x , d v , 2td t + xd x - ud u , 2td x -d u - 2xd v , 

4t 2 d t + 4txd x - 2(x + 2tu)d u - (x 2 + 2t)d v , e~ v {a x - au)d u - 4e~ v ad v ), 

is infinite-dimensional. Here a = a(t,x) is an arbitrary solution of the linear heat equation 
ott = ol xx . The elements of bi are the first prolongations of Lie symmetry operators of the 
potential Burgers equation Vt = v xx + v 2 that is equivalent to the linear heat equation under 
the transformation t = t, x = x, v = e v . Since this transformation is not point in the variables 
of the initial equation (u = ue v ), the structure of the space of potential conservation laws of 
the Burgers equation is more complicated than in case of the linear heat equation. Namely, the 
following theorem is proved in [17]. 

Theorem 1. The complete hierarchy of the potential conservation laws of the Burgers equation 
consists of one local conservation law 

Df,u — D x (u x + u 2 ) = 

and infinite number of the linearly independent potential conservation laws have the form 

D t (pe v )+D x ((p x -f3u)e v ) = 

parameterized by the linearly independent solutions (3 = f3(t,x) of the backward linear heat 
equation f3 t = Pxx- 

The potential systems of the second level of the Burgers equation (i.e., potential systems 
obtained with usage of potential conservation laws) are equivalent to potential systems of the 
linear heat equation. That is why, investigation of the second-level potential symmetries is 
reduced to investigation of potential symmetries of the linear heat equation. 
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6 Conclusion 



In the presented paper we construct potential symmetries and complete hierarchies of potential 
conservation laws of the Fokker-Planck and Burgers equations via reduction of them to the 
linear heat equation. A brief discussion is given for the potential symmetries of the linear heat 
equation obtained from the potential system associated with an arbitrary single conservation 
law of the linear heat equation. The case of characteristics 1 and x are studied exhaustively. 

Let us emphasize that the problem on potential symmetries of the linear heat equation have 
been solved only partially and is still open in the general statement. This problem includes 
investigation of potential systems ([8|) with arbitrary characteristics and of union of such systems 
with arbitrary families of linearly independent characteristics. In particular, it is necessary to 
find such conditions for families of characteristics that the associated potential system yields 
pure potential symmetries for the linear heat equation. 

The next important and interesting generalization of the obtained results is to realize the 
same program for arbitrary (1 + l)-dimensional linear parabolic equations. We strongly believe 
that the method used in [17] for construction of all possible potential conservation laws of the 
linear heat equation can be extended to this more general case. More precisely, it is known 
that the space of characteristics of the local conservation laws of a linear partial differential 
equation Lu = includes the functions of independent variables being solutions of the adjoint 
equation L*\ = 0. Our conjecture is that any local or potential conservation law of a (1 + 1)- 
dimensional linear parabolic equation is equivalent to one with the characteristic depending 
only on the time and space variables. The conjecture was already tested for the linear heat 
equation [17]. That is why, description of all possible potential symmetries of the linear heat 
equation will give a good hint to solve the similar problem for arbitrary linear parabolic (1 + 1)- 
dimensional equation. This will be the subject of our sequel paper. 
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